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0 Digital filter. 

® While wave digital filters (WDFs) can always be implemented by commercially available Integrated digital 
signal processors (DSPs), the efficiency of such implementations leaves room for improvement. Modified forms 
of WDFs are therefore derived that yield DSP implementations much more efficient than those previously 
available. The approach Is particularly ; suitable for lattice WDFs. The modified WDF-stmctures can be easily built 
in such a way that they possess all stability properties desired under signal quantization, thus not only absence 
of small- and large-scale limit cycles but even e.g. forced-response stability and related properties. 
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Digital Filter 



1. Introduction 

Wave digital filters (WDFs) /I/ are known to have many interesting properties /2/. The most important of 
these concern their excellent stability behavior, in particular with respect to alt aspects resulting from the 

5 nonlinarities that are caused by the signal quantizing operations required for carrying out 
rounding/truncation and overflow correction, and their good dynamic range performance. Other advanta- 
geous properties hold more specifically for voltage-wave digital filters (VWDFs). i.e., the type of WDFs 
usually considered and with which the present paper is also concerned: the small number of multipliers 
required, at least if the most appropriate structures are selected, and the small number of nonzero bits 

10 needed for implementing the multiplier coefficients. 

On the other hand, VWDFs require more adders than multipliers, at least if the most common type of 
implementation is adopted, i.e.. the one that Is the most economical in number of multipliers. This is not a 
disadvantage if a dedicated hardware is used, but can be undesirable If the implementation should be done 
by means of general-purpose digital signal processors of the types hitherto available. For these digital 

15 signal processors an individual computing cycle consists indeed always of one multiplication combined with 
an accumulation (addition). Hence, an individual addition requires a full cycle, i.e., a multiplication by 1 
followed by the addition itself. This can be a disadvantage for the implementation of WDFs, especially in 
the case of floating-point arithmetic. Recall tiiat In a WOF the arithmetic operators (adders and multipliers) 
are grouped in so-called adaptors. 

20 The purpose of the present method Is to show how the disadvantage just mentioned can be overcome 
in a wide class of WDFs. I.e., those for which only two-port adaptors are neeed. by reducing appreciably the 
number of cycles needed for carrying out the computations. This class includes the most attractive 
realization of lattice WDFs, i.e. those types of WDFs that have proved so far to be the most attractive ones 
in practice. It will be seen that the realizations thus obtained are such that they can make use to great 

25 advantage of "a feature available in several digital signal processors (see e.g. /3/). This feature allows one 
indeed to ensure that If an overflow occurs after a multiply-add operation the result is directly obtained 
according to a saturation characteristic. Due to the specific properties of WDFs tills guarantees not only 
suppression of overflow oscillations, but also forced-response stability and related properties /2.4-a/. 

Although we have so far stressed two-port adaptors. Uie method Is also applicable in the case of 

30 adaptors witii more than two ports. Considerable advantages are then still obtainable, but the reduction in 
number of cycles decreases with the number of ports of the adaptor. Note however that more than three 
ports are rarely needed in practice. 

The crux of the method consists In inserting, at selected locations of the original structure, appropriately 
chosen pairs of inverse multipliers and then to combine each such multiplier with the adaptor to which it Is 

35 adjacent. The pairs of Inverse multipliers do not necessarily conrespond to ideal transformers. Hence, we 
may no longer simply speak of port resistances or conductances (weights), but the two terminals of a same 
terminal-pair in the final structure may very well have different weights. This answers in a positive way a 
question that had been left open in a recent paper oh fully general passive and lossless digital filter 
structures /9/. It had indeed been mentioned there that it was "not known whether there exist situations in 

40 which further simplifications could result" by choosing distinct weights for the two terminals of a same 
terminal-pair. 

If not othenvise mentioned, terminology and notation will be as given in /2/. 



45 2. Reactances and all-pass structures Involving two-port adaptors 



2.1 Rrst-degree sections 

50 Since WDFs are based on using wave quantities rather than voltages and currents, realizing an 
Impedance or admittance is equivalent to realizing a reflectance; In particular realizing a reactance is 
equivalent to realizing an all-pass structur . A first-degree all-pass section is shown in Fig. 1a. Contrary to 
what we have priviously usually done /2/, we have indicated port conductances (weights) rather tiian port 
resistances, and we will do the same in corresponding later figures in this paper, such conductances 
(weights) being of course positive quantities. Consequently, for the multiplier coefficient 71 we have 
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71 = (G2-Qiy(G, + G2) , l7i| < I . (la.b) 

and the equations describing the two-port adaptor are 

bo = -7iao + (l + 7i)ai . (2a) 

b, = (l-7.)ao + 7iai . (2b) 
5 the subscript ! being given to y for reasons that will become clear later. We may assume 7*0 since 

othenwise the equations (2) become trivial and a problem of simplifying them does not arise. 

We now Insert pairs of inverse multipliers with coefficients kl , l/kl . kj , and l/k, as shown in Fig. 1b. 

Observe that the two multipliers of such a pair are placed directly in cascade. This does not change the 

signals ai and b, at the input port of the overall arrangement. We can combine the adaptor and its two 
70 adjacent multipliers into a new building block. Nj. as indicated by means of a broken line in Fig. 1b. This 

building block may be said to be frequency-independent just like an adaptor since it does not contain any 

delay. From Rg. lb we derive 

'5 r k'^a^ , b^ r • (3a;b) 

Hence. N| is described by 
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(fta) 



f 



(4b) 



where 



Yoo = -^1 ' ^01 = «i-Yl>»<l • (5a, b) 



35 We have two alternatives for making equal to unity one of the coefficients In each one of the equations 
(4a) and {4b). The first alternative is 



k^ = l/Cl^Yj^) , kj = ^' I' ^y^l = ^ -^1 ' (6a,b,c} 



which leads to 

7oo = -Yl . 7ol = 7lo i . 711 = 7/(1-7? ) • (7) 

The second alternative is 

45 

kj r l/(l*Yi) , s Yi/d-^Yi) » = (8a,b,c) 

50 

and leads to 

7oo = -71 . 7oJ = 711 = ' » 7Io = ('-7? )h\ ■ (9) 

Since |7i| < I. we can always ensure, if desired, that none of the coefficients In (4) is larger than I in 
modulus. This possibility is offered by (6) and (7) for |7i| ^ (V5-l)/2«0,6_18. and by (8) and (9) for I71I ^ (V5-I)- 
55 /2. In (6) and (8) we have supplemented the expressions for kl and kj simple expressions for kj /k| ; these 
show that we always have 
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Ik^/k'l ^ 1 . (10) 

There remain of course the two multipliers to the left of N| In Rgs. lb and 1c. These may obviously be 
combined into a single multiplier of coefficient kj /ki . In some cases this signle multiplier may even simply 
be dropped, but in others* this is not permitted. In any case, the implementation of N| requires only 2 
multiply-accumulate steps. A simplified overall representation of the structure of Fig. 1b Is shown in Rg. 1c. 

In addition to the port conductances Go and G|. the weights 



G* = k'^Gj^ , GJ = ^l^^i » (11a, b) 



75 are indicated in Fig. lb and 1c. Only G© corresponds to a true port conductance, while Gl and G| are the 
individual terminal weights of the two terminals of N| at its access on the left. These two terminals thus do 
not form a true port. The role played by Gl and G] will become clear in Section 3. 

It is obvious from what we have said that N| may be no longer termed a two-port» although it is a four- 
terminal (building) block or. simpler, a 4-pole block. 

20 

2,2 Sections of degree two and higher 

A situation in which the multipliers to the left of N| mentioned in Subsection 2.1 may not immediately be 
25 dropped is encountered in a second-degree all-pass section as shown in Rg. 2a. Note that we could there 
combine the two delays T/2 into a single delay T, but we have preferred keeping them separate for the 
sake of symmetry and for reasons as explained in /2/. We can again transfomi Rg. 2a in a way similar to 
what we have done previously, which leads to Rg. 2b. The port to the right in Rg. 2a is assumed to be 
transformed as in Rg. 1. with N| thus being as in Rgs. lb and 1c. The multipliers l/kl and k] , however, have 
2Q been shifted across the two delays T/2 and have been inserted into what is becoming Na. In addition to the 
weights G\ . and Gj given by (10) we now also have to consider terminal weights 
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40 



45 



50 



^2 ^2^^2 ' ^2 " ^2^^2 • ^^^^ 



A simplified overall representation of the structure of Fig. 2b is shown in Rg. 2c. 
The equations describing the left adaptor in Rg. 2a are 



= -^2^2 * (^*'''2^^3 ' (.13a) 
bj r (l.Y2>a2* ^2^3 ' (13b) 

where 

T2 = (G2-G,y{Q2 +G,) . |t2| < I. (14a.b) 

(14b) following again from (14a) because G}>0 and G2>0. We may assume 72^^0. Since (cf. Rg. 2b) 



®2 = ''l^Z ' ^2 ' ^1^2 ' (15a, b) 

55 83 = 63 s k^b^ , {15c,d> 



(1 3) gives rise to 
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b» = Y3,a' .^333*3 (16b) 



where 

10 



15 



= -Y^kJ/k^ , = (l*Y2)k*/k; . (17a, b) 

Y32 = (l-Y2)k^/k^ , = Y2*'2^^2 • (17c, d) 



These equations thus describe a new 4-pole block, N2, indicated in Fig, 2b by means of a broken line and 
shown In compact form in Fig. 2c. 

In (15). ki . and may be chosen arbitrarily, but ^22 cannot be influenced by the choice available. Like 
20 for N|, we have thus two choices for making equal to unity one of the coefficients in each one of the 
equations (16a) and (16b). i.e.. on the one hand. 



25 



k* r k*^/(l-..T2) , = (l-y^)kl , (i8a,b) 



which leads to 



Y22 ' -"^Z"! • "^23 = ^ (19a. b) 



^2 "1 



Y,,= 1 , Y = — • - . (19c. d) 

i-y; "^i 



33 y^^Z kV 

35 ^ ^2 
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(20a, b: 



and on the other. 

k^ = k*/{l*Y2) , k^ s ^2*^2 ' 
which leads to 



Y2 

These two choices still have to be combined with those available for N|, i.e., for k\ and k| . thus yielding 
altogether four choices. Due to I72I < I. taking into account (10). we can always ensure that none of the 
coefficients In (16) is larger than unity in modulus, and this independently of which one of the alternatives 
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for Nt has been adopted. The property mentioned is indeed achieved by (19) or (21) If 

i^^/d-Y^)! S IkJ/k* I or lY^/d-Yl)! ^ IkJ/k^l , (22) 
respectively. Note that from (18) and (20b) we also obtain 

kj/k^ = {l-Y^)kJ/k* and k^/k^ r , (23) 

respectively. These expressions together with (10) show that we have in all cases 

IkJ/k^l S 1 . (24) 

20 

The process described can be obviously extended to all-pass structures of arbitrary degree, say n (Rg. 
3a). We then proceed from the right to the left as we have done in Rg. 2. and we thus arrive at a structure 
of the form given In Rg. 3b. The derivation of Na will be exactly as that of Na, but with kj and kl taking the 
role of ki and k] etc. The equations describing the i-th two-port parallel adaptor, i - 1 to n. are given by 



10 



T5 



25 



30 



''2i-l = (i-^5«2i-2 * ^i'2i-l • (25b) 



where 

= (GrGH)/(Gi + G,,) < I. (26a.b) 
the Go to G„ being the consecutive port conductances. We introduce parameters kj to k^ and k] to k^ as well 
as new signal quantities a^ to a and b^ to b satisfying, for i = 1 to n. the equations 

®2i-2 ' ^i-l®2i-2 ' '*2i-2 ' **i-l*^2i-2 » (27a, b) 

^2i.l " '^1^21-1 ' ^2i.l = ^2i^l ' (27c, d) 

with 

45 

■•o ' ''S = 1 • «o = *o ' '»o = "o • (28a,b.c,d) 

^ This leads to 



55 



i-2 - ^2i-2,2i-2 ^2i-2 * ^2i -2 , 21-1^21-1 ' ^^^^^ 



^2i-l ' Y2i.i^2i-2 -^21-2 * ^21-1, 21-1^21-1 ' 
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where 

V2i.l,2i.2= ' ^2i.l,2i.I = • 

to ^ ^ 

These equations thus describe the 4-pole block N|. For the terminal weights of N| one finds 

G* = k[^tl^ , G? = kf G. , (31) 

at least for those at the left-hand side of N|. At the right-hand side of N| the terminal weights are the sam 
20 as the con-esponding ones at the left-hand side of Nj.,. They are thus those given by (31). but with I 
replaced by i-l. In view of (28a.b) this includes the case i = I. 

By appropriate choice of the k| and k] , we can make sure that, for each i = l to n, one of the coefficients 
in {29a) and one of the coefficients in (29b) becomes equal to unity. Each one of the N|, i = l to n, thus 
requires only two multiply-accumulate steps. The two solutions for which this holds are 



25 



30 which leads to 



''2i-2,2i-2 * ^2i-2,2i-l * ^ ' 



(33a, b) 



35 

and 

k[ = I<;.i/(l-^Yi) I kV c k*.^Y^/{l*Y.) , (34a, b) 

45 

which leads to 

Y2x-2.2i.2 » Y2i.2,2i-1 = ' ' ^"^'"^ 

2 „ 

^2i-1.2i-2 » ^ • ^ • ^2i-1.2i.l « ^ • f'^^"^' 
Furthermore, one of the two alternatives thus available for N| is such that it does not imply any 
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multiplier coefficient larger than I in modulus. Mor precisely, whil we have in general 2" possibilities for 
achieving the goal that none of the N| to N„ requires more than two multiply-accumulate steps, there is 
aJways one solution, and usually also only one, for which none of the final muitiplier coefficients is larger 
than one in modulus. This can t^e shown to be a consequence, among other things, of the fact that one has 
5 always either 



10 




or 



thus, in view of (26b) and (28a.b). 

75 

Uv /k'l i 1 for i s 1 to n. • (36) 



This is also the reason why it is never possible to arrive at a choice such as 721-2.21-2 = I- 
20 The two multipliers remaining at the left of Rg. 3b can be combined into a single multiplier m = !<„ /!<„ , 

for which we thus also have |m| ^ 1, with the inequality sign usually holding. In some cases this single 

multiplier may even simply be dropped, in others, it may in turn be combined in some fashion with another 

multiplier, thus leading also to a further saving. 

We conclude this section by listing a few useful relations that can be derived after some more or less 
25 tedious calculations. For the solution defined by (32) and (33) we obtain, using (26) and (31), 

^2i.2,2i-2 - 77771:7—; = 
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45 



SO 
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^2i-1.2i>I ^2i>2.2i>2 ,,0^ . v 

'2i-2,2i-2 ^2i-l,2i-l 

^[^i - "[.l^i.l ' ^'i - (1-^21-2. 2i-2 ^2i-l,2i-l>'^l.l- <^'^'''- 

Similarly, for the solution defined by (34) and (35) we obtain, using again (26) and (31). 



2 '^i^l 
^2i-2.2i.2 ■ ^2i-2,2i-2 ^2i-lt2i-2 " T7 



(40) 




s G[^^* , r ( 1 . ^2i-2>2i-2 j ^.^^ ^ (42, a, b) 

^21-1,21-2 



3. Realization of circuits with more than 2 terminals 
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WDFs in lattice configuration hav so far turned out to be th nnost attractive ones from a practical point 
of view 12J, A lattice WDF can be built by means of two WDF all-pass structures. For these, the results of 
Section 2 are immediately applicable. In fact, one can make use of these results In different fashions since 
any all-pass function can be realized by means either of a chain connection of unit elements or a cascade 

5 of individual sections of degree, say. one and two. in the former case, we obtain a WDF realization in the 
form of Rg. 3a, thus leading to the structure of Rg. 3b. In the latter, we are lead to making use of th 
results described In relation with Rgs. 1 and 2. In all cases, however, there obviously will be. for each 
complete all-pass structure, only one multiplier that has to be implemented in addition to those inside of 
building blocks such as N, to Nn in Figs. 1 to 3. i.e.. in addition to those implemented in form of muliply- 

;o accumulate operations. 

A configuration of a lattice WDF for the case that only one input and one output terminal is used is 
shown in Rg. 4a. It comprises two branches of all-pass transfer functions Si and S2. while A, and B2 are th 
(steady-state or z-transform) Input and output signals, respectively. (The factor 2 at the output is irrelevant; 
it has been Included in order to be in confomnity with Rg. 23 In 121) Let Sl and Si be those functions 

75 corresponding to Si and S2, respectively, that are obtained by the method explained in Section 2 if we 
ignore the remaining multipliers mentioned In the previous paragraph. 

Clearly, the structure of Fig. 4a is equivalent to that of Fig. 4b where mj and ma are the coefficients of 
the multipliers just referred to. Clearly, the structure of Fig. 4b is equivalent to those of Rg. 4c and 4d. In 
these, the multipliers m, and ma, respectively, may again usually be dropped, and the remaining multiplier 

20 may be combined with the adder into another multiply-accumulate operation. Note that the discussion given 
here does not consider scaling requirements. 

The situation is. obviously, very similar if the full two-port lattice WDF (Fig. 5a) is to be implemented. 
One of the two main possibilities then available is shown In Rg. 5b. 

There are also other WDF structures for which the method is directly applicable. This Is the case in 

23 particular for a WDF (Fig. 6b) obtained from a reference filter in form of a cascade of unit elements (Fig. 
6a). We may in this case again start from the end with the port conductance Go, similar to what wejiave 
done in Section 2, and then proceed to the other end, which leads us to the structure of Fig, 6c. Note that 
in Rg. 6 we have adopted a numbering of the G,. i=0 to n + l, which is different from the one we have 
usually chosen in the case of doubly-terminated filters /2/, but which is more in conformity with the on 

30 used in Figs. 1 to 3. 

The process is precisely as that described in Section 2, except at the beginning. Indeed, assuming that 
a multiplier in cascade with an externally accessible terminal is inrelevanti we may introduce multipliers l/ko 
and ko in cascade with the input and the output terrnihal. respectively, of port 0. If we choose ko =ko =1 (cf. 
(28a.b)) everything will be as in Section 2, but if we accept ko *ko we have one additional degree of 
35 freedom. This may be used to make one further coefficient (e.g. in N,) equal to unity. 

The methods explained in this section may. of course, also be of interest in the case of circuits having 
more than one input and more than one output terminal. 



40 4. Stability considerations 



4.1 Ideal lossless building blocks 

45 It is known that the most important advantages of WDFs are their excellent stability properties, 
especially those under the various nonlinear conditions resulting from the unavoidable quantization, oper- 
ations needed for the signal quantities. Due to the way we have derived the new circuits described in 
Sections 2 and 3, from conventional WDFs it follows that the former must inherit all stability properties of 
tfie latter. We will briefly examine the mechanism behind this observation. We do this first under th 

50 assumption of ideal lossless building blocks N,. The case of nonideal blocks N, will be considered in 
Subsection 4*2. 

Consider thus a 4-pole N,, I €{1,2 n}, of Rg. 3. It Is described by the equations (26) to (31), the latter 

to be used, for a given N|, also witii i replaced by i-l. In order to simplify the writing, we consider specifically 
the case i=2. in which case the equations just mentioned may be replaced by (11). (12), and (14) to (17). 
55 (Note that i = I would not be appropriate if one wants to be general, since for N| the specific conditions (28) 
hold.) 

The instantaneous power (pseudopower) absorbed by the original adaptor with coefficient 72 is given by 
P2 = (G,^ + Gaai ) - (G|b| + Gab? ) (43) 
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It is known (and can be verified directly) that, due to losslessness. we have p2=0. Using (11), (12). and 
(15). we obtain from (43), 

P2 = ^^1^2^* ^2^3^' - ^^l''2^* ^2*'3^' ' ^^^^ 

We note that in (44) each one of the four signals a2 . 83 , bz , and bs appears squared and multiplied by the 
weight of the particular terminal to which the signal refers. Under ideal linear conditions, i.e.. if (16) holds. 
(44) yields again p2 =0. assuming of course that (11). (12), (14a), and (17) are satisfied. 

Mathematically speaking, the main difference between (43) and (44) consists in the fact that while In 
(43) the second paranthesis comprises the same weights (i.e., G| and G2) as the first one, this is not the 
case for (44). It is known however /2.9,10/ that this does not affect any of the stability proofs that have been 
provided for WDFs. This implies that the same type of quantization rules that guarantee, If applied to the bj. 
a specific type of stability In the original WOF. will guarantee the same type of stability in the modified WDF 
if applied to the b\ , with I ranging over an appropriate set of integers. 

In order to be more specific, let us designate by b jq the value resulting from bl by quantization. There 
will be no observable small-scale or large-scale limit cycle if quantization is carried out In such a way that' 
|b [ql ^ |bi I for all relevant values of i. If the circuit is such that it cannot sustain unobservable periodic 
oscillations under ideal linear conditions and if the rule adopted for overflow corrections is such that a 
simple sign inversion is excluded, there will be no unobservable limit cycle either. The conditions 
mentioned allow us also to guarantee stability under looped conditions and to state equally that the limit 
cycles superposed to the output signal if the input Is an arbitrary periodic signal will be very small. 

Of particular interest finally is the fact that forced-response stability and related properties /4-8/ can be 
guaranteed by simply requiring that for all relevant i*s the value of b [q Is obtained from the corresponding 
b'i by adopting for overflow correction e.g. either simple saturation or a triangular overflow characteristic with 
slopes of i45 * . The first of these possibilities is directiy available in several digital signal processor (e.g. in 
the TMS 320 /3/) at tiie completion of each multiply-accumulate step. Thus, since the computation of any 
b [q requires just one such step, using the option available in such a digital signal processor automatically 
guarantees forced-response stability and related properties. 

It should be stressed that all this holds for any type of signal representation, in particular thus not only 
for fixed-point arithmetic but also for floating-point arithmetic. This is particulariy important since any 
floating-point digital filter that cannot be built strictly without parasitic oscillation can always sustain a 
parasitic oscillation involving the highest possible value of the exponent /11,12/, In practice thus a parasitic 
oscillation of high amplitude. 

It should finally be recalled that an expression such as that appearing in tiie right-hand side of (44) 
corresponds indeed to the definition of the power absorbed In the case of a class of digital filters that 
appears to be the only one that offers all the same good features as conventional WDFs. yet is in a sense 
somewhat more general /9/. It had been shown in /9/ that this extended class can always be obtained from 
conventional WDFs by precisely tiie type of transformation described in Section 2. Since this type of 
transformation Is of rahter trivial nature, the digital tilters described in this paper should still appropriately be 
called "wave digital filters". 

In /9/ the question had been left open whether the transformation just referred to can lead to structures 
offering true advantages over those obtainable by the more conventional WDF approach. The results of the 
present paper show that Uiis question can now be answered in the affirmative. 

4.2 Lossy building blocks 

In a conventional WDF an adaptor can be easily made strictiy lossless, at least under the assumption 
that all additions and multiplications are carried out exactixy. I.e., that the circuit is strictly linear. This is due 
to the fact that e.g. for a two-port adaptor described by (13) strict losslessness is fulfilled for any value of 72 
provided we adopt port weights G| and G2 that satisfy (14). Since the latter requirement is trivial to meet, we 
can. in particular, adopt for yz any value expressible in a binary representation witii finite number of bits. In 
this case, all four coefficients in (13). In particular thus also the coefficients (I-72) and (I + 72). are 
expressible in the ssam way. 

This same simple property does not hold for the blocks N|. 1 = 1 to n. say. in an all-pass section of the 
type of Rg. 3b. In order to show this we simply consider N|. Since Go =Go =Go we obtain, for i = l. from 
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(37) and (40) 

- lii = 1 and . y y . = i , (45a, b) 

00,^ oo oolo 

•11 

respectively. Thus 700 and 711 (first equation) as weil as 700 and 710 (second equation) have to satisfy an 
expression that usually cannot be fulfilled if these coefficients have to be expressed in a finite number of 
binary digits. 

In order to solve this dilemma we observe that for all stability aspects discussed in Subsection 4.1 to 
hold it Is merely required that the blocks Ni are passive, i.e.. such that pi S 0, 1 = 1 to n. and this for all values 
of the input signals of N,. pi being defined as given by (44) (for j = 2); strict losslessness. I.e.. pi = 0 Instead 
of pi ^0. is not needed. 

We first consider N,, more specifically the choice leading to (7). Let 7ooq and 7iiq be the quantized 
values of 700 and 7h. respectively. We may assume the value of Go to be given, but Gl and G may have to 
be replaced by new values, which we designate by Q ;q and G Vq, respectively; suitable choices for 
these still have to determined. Let piq be the power absorbed by N, under these new conditions. We have 

Using (4) (but with 700 and 711 replaced by 7ooq and 7Pq. respectively) and 701 = Yio=i. this yields 



36 



40 



45 



SO 



Any choice of 7ooq. 7uq. G and G yq for which piq i 0 for all ao and al is acceptable. 

We will now show that the problem thus formulated Is Indeed solvable. For this, let us choose G Jq 
and G Vq ^^^^ 

g' - (l-Y Y , )G , GV = - G„ , (aaa,b) , 

^Iq - ^"^^ooq ^llq^^'o * ''Iq 7 ° 

llq 

these expressions conresponding in a sense to (39b) and (38b), respectively, with i = l. In order to lead to the 
desired stability results, these expressions must yield positive values of the weights G {q and Q 'jq. For 
this, observe that In view of (lb) and (7) we have 700711 < 0. and we may thus assume the quantization to be 
such that the signs of 7 00 and 7n are not reversed. I.e.. in particular, that 
Tooq 7iiq ^ 0 , -7o<H/7iiq ^ 0 . (49) Substitution of (48) in (47) yields 



' ^ooq > «2 

r 

llq 



223 )G ' (50) 

o 0 



and thus requires 

7 §oq - (7ooq/7.1q) ^ I . (51) 

which, in view of (49). is equivalent to 

7 Soq + l7W7IIq| ^ i ■ (52) 

Comparison with (45a) shows that It Is easy to carry out a quantization of 7 00 and 711 such that (61) is 
fulfilled. All requirements are then m t if G [q and G Vq are chosen according to (48). 
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Next we consider again N|, but with the choice leading to (9). From (46) we now obtain 



- 2^Wo*^loqSq>S<- 

10 

By an approach similar to that used before one finds that pic^ 0 for all ao and aj if e.g. we quantize and 

7io in such a way that 

y ooq - 7ooq7toq ^ L (54) 

and if we select G Jq and G *;<, such that (cf. (41b) and (42a) for i = 1) 

IS 

GV = -(y /Y, )G . G! = GV ♦ . (55a, b) 

Iq ^^ooa 'loq' o- Iq Iq o 

Assume next that the circuit involves also N2. Having performed a quantization of the coefficients of N| 
and having thus also chosen G Jq and Q 7q, we may now approach Nz. We assume first that for Nz the 
choice leading to (19) has been adopted. Let 722q and Tsaq be the quantized versions of 722 and 733. 
respectively, let Q ^q and G 'i^te corresponding choices for G2 and G2 * respectively, and let p2q be the 
value thus replacing p2 given by (44). We have 

26 

P2q ' "^iq -2'* "^2, " '^Iq ^'z^ " ' 

30 

From this we obtain, using (16) (with 722 and 733 replaced by 722q and 733q. respectively) together with 
(19b) and (19c). 



35 



40 



45 



''2q = (^iq- >22q^'lq-'^2q'-2' * < '^2q-'^lq-i3q<^2q ^ ^ 3^ 

- 2<^22q'=lq*^33q'^2q'«2'3 ' <57) 

where G jq and G 'fq may be assumed to be known from quantizing the coefficients of N|. Any choice of 
722q> Y33q. G 2q> ^ 24 which paq^ 0 for all 32 and as is acceptable. 
One such choice is obtained if we put 

^2q = «^-^22q^33q>Clq ' ' " 22q^ 33q ^'^iq ' (^8«''>' 



50 these expressions corresponding, in a sense to (39b) and (38b). respectively, with i = 2. This leads again to 
positive weights G ^ 2q> indeed it follows from (14b). (19a), and (19d) that 722733 < 0 so that we 

may assume 

Y22q 733q < 0 . 722q^733q < 0. (59) 

Inserting (58) into (57) leads to 
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from which we derive the requirement 

(61) 



Hence, a quantization according to (61 ) together with a choice of G gq ^ aq according to (58) meets 
all the requirements. 

If for N2 we adopt the choice leading to (21), it turns out that (57) Is replaced by 

10 

P2P = ^^U ' 'Iz.'l'^ • ^W5c,)^2' - ^'^2. - 'U'^'^" 

We now have to perform the quantization of 722 and 732 as well as the choice of G gq and G such a 
^ way that paq^O for ail a2 and 33 . A possible simple solution for this can be found in a similar way as before, 
i.e., by combining the requirement 



25 



>22q -^22q ^32q ' ^Iq^^lc 



(63> 



with the choice 

r.. .'t^ /V ^^• g! = C' + . (64a. b) 



t^Sq = -^^22q/^32q>'^lq ^2q= ^Iq * ^'i^ 



If n >2, the process can be continued until we reach N„. Thus, for any is {1.2 n} and with G and 

3^ G V.i.q known, we must quantize the coefficients of N| in such a way that piq S 0 for all choices of a 21-1 
and a 21^. with 



40 



A possible solution for this is immediately obtained as generalization of our results given above. Thus, if for 
Ni we adopt the choice leading to (33) we may first quantize 72h2.2w and 72h^h according to 

^ 2 ^2i>2.2i-2,q ^ ^i-l>q (55) 

Y2i-2,2i-2,q * ^ , 

so and then choose G |q and G 7q according to 

' ^l-^2i-2.2i-2.q^2i.l.2i-l.q^=i.l.q ' ^"^^ 

55 

= -C^2i.2.2i.2.q/^2i-1.2a.l.q>'^i.l.q * ^""^ 
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Similarly, if for Nj we adopt the choice leading to (35) we may quantize yai-a^ and tzh^w according to 

G" , 

2 Y Y S (67) 

^2i-2,2i-2,q ^ ^2i-2 , 2i-2 , q '2x-l , 2i-2 , q , 



and then choose G \q and G *j*q according to 



^'U - -<^2i.2.2i.2.q/^2i-1.2i.2.q>'^i*.l.q ' 
^Iq - «^I.l.q * ^iq ' 



if in such a process an optimum solution (say. with respect to the transmission properties of the circuit) 
is desired for the quantized coefficients, one should of course repeat the full cycle from I to n as often as 
needed. Note that we could then encounter, for one or more of the ie {1,2,...,n}, a case for which one of the 
two multiplier coefficients may be chosen equal to zero. We then also choose the other coefficient equal to 
zero and we restrict this situation to the case leading to (33), thus to 

72i-2.2l-2,q ~ 72i-i.2W.q = 0 

and 



«^iq - ^i-l.q • "^iq ' «^i-l.q ' • 

The problem of determing Ni thus then becomes trivial. 

After having completed the quantization of the coefficients of the blocks N|. i 1 to n, there may remain 
the problem of choosing l/kn and k'n (cf. Rg. 3b), or simply that of choosing 



in s k"/k' 
n n 

A need for this does not arise if m may simply be dropped. In other cases such as in Fig. 4c the actual 
problem is that of finding a quantized value for ma/mi, where mi and ma in turn may be products of factors 
such as m, but this is not a question of stability. 

There may be cases however where passivity from input to output must be guaranteed. In such cases it 
is best to consider the corresponding transfer function. Let thus A 2n4 and B be the steady-state (or z- 
transform) quantities corresponding to a 2rH and b in-i* respectively (Rg. 3b). If all blocks N| to N„ are 
designed as discussed above, i.e., if these are all passive, we have /9,10/ 



i.e.. for the transfer function S. 



IS!2 S Ig;/G;;! -here S = B^„,i/A^„.i . 
Note that in our stability analysis we have specifically examined the cas of the structures of Figs. 1 to 
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3. it is obvious that similar considerations hold for other structures, e.g. for the one of Rg. 6. 

In the above procedure of quantizing the multipli r coefficients in such a way that the blocks N|. I-I to n, 
remain passive, we have, for each case considered, given a specific solution by which the desired goal can 
be achi ved. It should be str ssed that other solutions may very well be feasible. Thus, considering the 
general case of Ni, adopting the choice leading to (33). we may replace (66a) by 

«^iq = -<Y2i-1.2i.l.q/^2i.2.2i-2.q>'5;L-l.p ' ^^'^ 

It can be shown indeed that perfonnning the quantization according to (65) and selecting G {q and G jq 
according to (69) and (66b). respectively, also meets all the requirements. 

However, such an alternative cannot be better than the one given by (65) and (66). Indeed, using (65) 
with i replaced by i + l. we see that the quantity K, = G ;VG Jq should be as large as possible in order to 
facilitate as much as possible the quantization of N^i. One verifies that the ratio of the values of Ki obtained 
by choosing G jq according to (69) and (66a). respectively, is equal to 



20 . ^2i^2,2i-2>q j ^^1^ 

^Y2i.2,2i-2,q Y G" 

2i-l,2i-l,q i-l,q 

which according to (65) is ^I. 
25 The situation is similar if for Hi the choice leading to (35) is considered. An alternative solution is then 
obtained if we replace (68b) by 

' -n'-l,q/^^2i.2.2i-2.q^2i.l,21-2,q) ' 

but again, it is never better than the one given originally. The reason for this is the same as that given in the 
previous paragraph. 



5. WDFs involving adaptors with more than two ports 

In certain types of WDFs one has to make use of adaptors with more than two ports. The question thus 
arises whether the present approach can be extended to such cases. 

In principle this question can be answered in the affimnative. It is Indeed always possible to introduce 
pairs of inverse multipliers as we have done In Section 2 and to make use of the freedom thus gained in 
order to introduce simplifications for the realization of the adaptors. Hence, one can take advantage of this 
possibility also in the case of adaptors with more than two ports. 

In general, however, one cannot expect to arrive at as substantial a gain as in the case of two-port 
adaptors. The reason for this is that the number of entries of a matrix describing an adaptor increases with 
the square of the number of ports, while the number of degrees of freedom gained for an adaptor increases 
only linearly with the number of its ports. 

It is thus clear that beyond the two-port adaptors the method appears to be most attractive in the case 
of three-port adaptors, which are also the ones that are most important in practice. We will not examine this 
in detail in this paper, but simply point out a few aspects. 

We will have to expect that in the case of a three-port adaptor the resulting 6-poIe block will be such 
that 2 multiply-accumulate steps will be necessary for each one of the 3 output signals (except for the 
output signal at a reflection-free port, which needs only one multiply-accumulate step). A difficulty arises in 
this case with respect to implementing a simple saturation characteristic, at least if the only possibility of 
doing it is to apply it after each multiply-accumulate step individually and if the three values to be added do 
not all have the same sign. A way of avoiding this difficulty is offered if none of the coefficients involved is 
larger than unity in modulus and if there exists a possibility to an^ange the order of the computations in such 
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a way that the first accumulation step concerns two numbers of opposite sign. Even without the latter type 
of facility, however, overflow stability (as apposed to the more stringent forced-response stability) can be 
shown to remain guaranteed. 

5 
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35 Figure captions 

Fig- 1 (a) A Ist-degree all-pass section 

(b) A section derived from (a) by Inserting pairs of inverse multipliers at appropriate locations 

(c) Simplified representation of (b). 

40 Rg. 2 (a) A second-degree alt-pass section 

(b) A section derived from (a) by inserting pairs of inverse multipliers at appropriate locations 

(c) Simplified representation of (b). 

Rg. 3 (a) An all-pass structure of degree n 
(b) A structure derived from (a) by inserting, at appropriate locations, pairs of inverse multipliers and 
45 combining them (except for the two multipliers at the left) with the adaptor to which they are closest 
Fig. 4 (a) A lattice WDF with one input and one output terminal 

(b) Corresponding structure obtained by realizing the all-pass functions S| and Sz as explained in Section 2 

(c) and (d) Equivalent structures derived firom (b). 

Rg. 5 (a) Full two-port lattice WDF 
50 (b) Example of a structure derived from (a) by realizing the all-pass functions S| and Sz as explained in 
Section 2. 

Rg. 6 (a) A reference filter in form of a cascade of unit elements 

(b) A corresponding WDF 

(c) A structure derived from (b) by means of a method analogous to that described in Section 2. 
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Claims 

1 . Digital filter obtained by modifying a wave digital filter (WDF) characterized in that the original WDF 
is modified by inserting, into appropriately selected interconn cting leads, appropriately chosen pairs of 

5 (mutually) inverse multipliers (kii l/k|...) and by incorporating, wherever possible, each such multiplier (ki; 
l/k|...) into the adaptor to which it is most adjacent. I.e., to which it is either directly adjacent or from which It 
is separated at most by a delay element wherein a structure resulting from an adaptor by Incorporating into 
it the multipliers (ki; l/k|...) to which it is most adjacent being referred to hereafter as a modified adaptor. 

2. Digital filter as claimed in claim 1 characterized in that the multipliers (kt; l/k,...) are chosen in such a 
10 way that the number of multiply-add operations required by the modified adaptors are as small as possible. 

3. Digital filter as claimed in claim 2 characterized in that for any modified adaptor resulting, from a two- 
port adaptor the number of multiply-add operations required is exactly 2. 

4. Digital filter as claimed in claim 3 characterized in that for any modified adaptor the values of the 
multiplier coefficients involved are both at most equal to unity. 

15 5. Digital filter as claimed in one of the claims 1 to 5 characterized in that the modified WDF Is a lattice 
WDF, 

6. Digital filter as claimed in one of the claims 1 to 5 characterized in that the modified WDF Is such 
that it is obtained from a WDF whose adaptors are all two-port adaptors. 

7. Digital filter as claimed in one of the preceding claims characterized In that the arithmetic operations 
20 in the modified adaptors are carried out in such a way that the output signals are detrermined according to 

a saturation characteristic. 

8. Digital filter as claimed in one of the preceding claims characterized in that in the modified adaptors 
Involved, the multipliers (kt; l/k|...) are quantized In such a way that these modified adaptors remain passive 
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